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X New address: Institut fiir Quantcnoptik, Leibniz Universitat Hannover, 30167 Hannover, Germany. 
§ New address: Department of Physics, Purdue University, West Lafayette, IN 47907, USA 



Resonant five-body recombination in an ultracold gas of bosonic atoms 



2 




Figure 1. TV-body scenario in the region of negative two-body scattering length a. The 
lower panel shows the iV-body binding energies as functions of the inverse scattering 
length. = (fiK,) 2 /m is the trimer binding energy for resonant interaction. The 
dotted, dashed and solid lines refer to the ground states of three-, four- and five-body 
cluster, respectively. Excited states are not drawn here. The pictorial representations 
of the three-, four-, and five-body bound state illustrate one of the possible spatial 
arrangements of atoms, given that an iV-body cluster state has N(N — 1)(N — 2)/6 
possible combinations of the three-body sub-clusters. The upper panel shows the 
relative losses from the atomic sample for a typical set of measurements, performed 
with a hold time of 100 ms for a sample of 10 5 atoms at 100 nK in a trap with mean 
trap frequency of about 30 Hz. 



1. Introduction 

Few-body physics with ultracold atoms has emerged as a new research field combining 
concepts from atomic, nuclear and condensed-matter physics. A growing number of 
experimental and theoretical studies have been focused on both the fundamentals of 
few-body phenomena [1] and the connections with many-body systems j2j [3]. The 
cornerstone of recent experimental advances is the control of interactions in an ultracold 
atomic gas, offered by magnetically tuned Feshbach resonances jl]. In particular, 
the tunable s-wave scattering length a allows to access the regime of resonant two- 
body interactions. Here, the system is governed by universal behavior, independent 
of the short-range details of the interaction potential. The paradigm of universality 
is Efimov's solution to the problem of three resonantly interacting particles [5]. Once 
the intimate connection between Efimov states and three-body recombination has been 
established [IJ|6j[7], resonant loss features became the fingerprint of Efimov physics in 
experimental studies with ultracold atoms [8l-f2~0]. 

Advances in three-body physics led to intriguing questions on the generalization of 



Resonant five-body recombination in an ultracold gas of bosonic atoms 



3 



Efimov's scenario to more particles and on the existence and observability of universal 
few-body cluster states [2TH26] . iV-body cluster states known as Brunnian states 
exist in a range of interaction where no (N — l)-body weakly bound subsystems are 
present [27J|28]. The general connection, however, between cluster states and the three- 
body Efimov states has remained an open issue. It was soon realized that no "true" 
Efimov state jjj] with N > 3 exist, because of the quite different scaling and threshold 
properties of the cluster states [29]. However, other approaches to extend universal 
Efimov theory to larger systems have been pursued along different lines [30H32] . The 
development of accurate descriptions of four-boson systems [33H35] demonstrates the 
existence of four-body states tied to each three-body Efimov state. In a major extension 
of Efimov physics, Ref. [36] predicts the existence of a family of cluster states tied to 
an Efimov trimer. 

According to Refs. [34T - l3~T] . the binding energy of the cluster states follows universal 
scaling laws, which are directly connected to the Efimov effect. Figured] shows the 
calculated energy spectrum for the ground states of three-, four- and five-body clusters 
(lower panel). The corresponding experimental observables are loss peaks in the atom 
number (upper panel), which appear at values of the scattering length a_, _, and 0,5 _, 
where the three-, four- and five-body states cross the free atom threshold, respectively. 
The resonant values of the scattering length are predicted to be universally connected 
by the simple relations a^- = 0.44(1) a_ and as ; _ = 0.65(1) a^- [37]. A more recent 
study [38] has theoretically explored the four-boson resonance and has determined 
its position with greater accuracy. For four-body states, the universal relation has 
been confirmed in experiments [12l[l9j[20j[39] and the four-body recombination rate 
has been measured [39] and calculated within the hyperspherical framework [34]. A 
straightforward extension to five-body systems is not currently possible for experiment 
nor theory. The experimental challenge is to discriminate the five-body recombination 
signal against a strong background resulting from fewer-body processes. The numerical 
difficulty of the scattering few-body problem grows exponentially with the number of 
particles making the description of five and larger systems beyond current theoretical 
capabilities. 

This article presents a combined theoretical and experimental study of universal 
few-body physics up to five-body states. We present strong evidence for the existence 
of an Efimov-related cluster state of five identical bosons and we provide quantitative 
results for the corresponding five-body recombination rate. Our results highlight a new 
level of understanding concerning few-body physics and its experimental manifestations 
in ultracold atomic quantum gases. 

|| The "true" Efimov effect refers to the appearance of an infinite number of iV-body bound states, 
which have a discrete scale invariance and which exhibit well-defined thresholds given by the (N — 1 )- 
body subsystem. 
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2. Theoretical approach 

The theoretical analysis of iV-body recombination processes requires the description of 
the iV-body scattering continuum. The hyperspherical framework has been successfully 
applied to describe recombination processes for N > 3 (341140] . In this framework, the 
Hamiltonian is diagonalized adiabatically as a function of the hyperradius R, which 
describes the overall size of the system, leading to a set of coupled one-dimensional 
Schrodinger equations. At ultralow temperatures and large scattering lengths, iV-body 
recombination events are mainly controlled by scattering processes with incoming flux 
in the lowest iV-body scattering channel and outgoing flux in deeper loss channels. 
The coupling to the deep channels is assumed to remain approximately unaffected 
as the scattering length and the collision energy are varied throughout our regime of 
interest. Thus, all the relevant information comes from the analysis of the lowest TV-body 
potential curve corresponding to the incoming scattering channel. 

However, the extraction of the hyperspherical potential curves and couplings 
becomes computationally unfeasible as the number of particles increases. Currently, the 
numerical technology allows for the calculation of potential curves of three and four- 
body systems but no tractable method has yet been implemented that can calculate 
both the potential curves and couplings for N > 4. However, the trapped energy 
eigenvalue spectrum of a five-body system can be accurately obtained with current 
technology. Thus, the present study extracts the relevant recombination information 
from an analysis of the trapped spectrum. 

Our starting point is the hyperspherical description of the ultracold iV-body 
recombination rate 140 1, 



\k 2 J 



where fM^ — m/ yN with m being the atomic mass, k = {2^ N E/hfy/ 2 is the incoming 
hyperradial scattering wavenumber and Sqq is the diagonal element of the S-matrix 
for the lowest channel (00) in the J n = + symmetry. For purely elastic scattering 
I^qq | 2 = 1 and the recombination rate is zero. In the limit in which every iV-body 
collision leads to losses, \Sqq \ 2 = 0. Taking the thermal average in the full loss 
temperature T, one obtains the unitary limit for TV-body recombination at low energy: 

. rp / % 2 x3(JV-l)/2 

{Ln)t = ^ (_y . (2, 

Based on potential curves computed for the three- and four-body cases, we expect 
that the five-body potential curve should have the topology depicted as a solid curve in 
Fig. [21 i.e., the lowest potential curve exhibits a barrier that separates the inner region 
(small R) from the asymptotic scattering region at large R. Note that the effective 
mock-centrifugal barrier j?T] in the lowest iV-body continuum channel at large R is 
guaranteed to have the form (for finite a), U(R) ->■ h 2 (3N - 4)(3iV - Q)/2/i N R 2 . For 
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Figure 2. Schematic representation of lowest five-body hypcrsphcrical potential curve. 
The solid curve represents the free-space hyperspherical potential while dashed curves 
represent hyperspherical potential curves in the presence of a harmonic potential. Solid 
lines represent the inner and outer WKB phases. Inset: Five-body trapped energy 
spectrum in the region where the bound five-body state crosses the lowest trapped 
states. 



this potential curve topology, an extension of the semiclassical (WKBJ) treatment of 
Berry [12] to include the decay to the lowest channels yields 



/, , „ ,o x e 27 sinh(2?7_) , . , . . , 

(1 " ISool 2 ) = — — j- I > A(77_, 7 ,0in , (3) 

l cos^ 0; n + smn J]_ 

where 0j n is the WKBJ phase for the inner allowed region, 7 is the WKBJ tunneling 
integral in the barrier region, r\_ describes the decay to deeper, non- universal, channels 
and is treated as a fitting parameter, and A(r)-, 7, 4>- in ) ensures the proper normalization. 
Thus, the determination of 7 and 0; n for the relevant range of energy and scattering 
length gives an approximate description of five-body recombination. 

The novelty of our approach is the determination of 7 and <f)- m from an analysis 
of the trapped spectrum. The five-body trapped energy spectrum is obtained using 
a correlated Gaussian basis set expansion [57] ; for details see Appendix A In the 
region where the five-body resonance occurs, the spectrum exhibits a series of avoided 
crossings between the five-body states bound in the inner region and the outer trap 
region states (see Fig.[2]). In the WKBJ approximation, the quantization condition for 
the trapped potential curve (dashed curve in Figj2]) is ^ tan((/> in ) tan(0 out ) = 1, where 
A = e _27 /2. For collision energies below the barrier local maximum and away from the 
avoided crossings, the allowed energy eigenvalues occur when (f> a (E,a) ~ ir(i + 1/2), 
where i is an integer and a = in, out. The phase 0j n in the four-boson case near 
the resonance energy is known from our previous work [20] to be well- described by 
(pin ~ 0in,o + H a A*vdw) + caE / (r vdW E vdW ) , where 4>m,o, b and c are fitting parameters 
and r v( jw and E vd w are the van der Waals radius and energy, respectively, as defined 
in [1]. Using this simple form in the five-body case and imposing the eigenstate condition 
for bound states, the values of (pm.o, b and c for = 5 are extracted. 
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Next, an analysis of the spectrum at the avoided crossings (see e.g. Fig.|2]) 
determines 7, as explained in detail in Appendix B The relevant avoided crossings 
occur when in »s out ~ 7r(z + 1/2). For narrow avoided crossings (A <C 1), the 
quantization condition reduces to S(j) in 5(j) out ~ ^ where <p a = n(i + 1/2) + S(p a . Right 
at the avoided crossing, the energy difference between the two states (2AE) is related 
to A, namely as 5<f) a ~ (d<f) a /dE)AE. The quantization condition thus reduces to 
A ~ 2AE 2 (d(f>i n /dE)(d(f) ou t/dE). Consequently knowledge of </> in , </> out and the energy 
avoided crossings allows the tunneling 7 to be determined. At large R, interactions can 
be treated perturbatively leading to an hyperspherical potential curve valid at large R 
that determines </> out . 

By changing the trapping confinement we can change <j) out and by introducing short 
range three-body forces we can modify <p- m without affecting the barrier. This allow us 
to explore how 7 depends on both E and a. Interestingly, at low energies our numerical 
results are nicely fitted by the formula: e~ 27 oc (E/E vdW ) 5 \a/r v< iw\ 9S , which is in good 
agreement with the predicted threshold behavior (1 — | -S'ocT 1 2 ) ^ -E^H 10 jSQ] - 



3. Experiment 

We prepare an optically trapped sample of cesium atoms in the lowest sub-level of the 
electronic ground state under similar conditions as described in Ref. [13]. The final 
evaporation process, performed at a magnetic field of 894 G (a = +285 a , where a 
is the Bohr radius), is stopped before the onset of Bose-Einstein condensation. The 
sample is then adiabatically recompressed to avoid further evaporation from the trap. 
At this point, the sample contains about 6 x 10 4 atoms at a temperature T = 78(3) nK 
and the confining optical potential has a mean trap frequency u = 2n x 36.2(2) Hz, 
which results in a peak number density of 4.2 fj,m~ 3 and a peak phase space density 
close to unity. The wide s-wave Feshbach resonance with its pole at 786 G offers ideal 
tuning properties [331 HI], superior to the low-field region investigated in our previous 
works 0[39]. 

We measured the decay of the atom number in the region of negative scattering 
length (from — 500 ao at 863 G to — 200 a at 873 G), where the four- and five-body 
recombination resonances are expected. After the recompression stage, we tune the 
scattering length to its target value, and we measure the atom number after a variable 
hold time by absorption imaging. After 100 ms, the typical loss fraction is around 10% 
at — 300 ao and almost 35% at — 450 ao- 

The time evolution of the number M of trapped atoms and temperature T are 
determined by the different A^-body loss processes and can be expressed [35] in terms 

of a system of coupled differential equations, 

+00 

M/M=-Y,L N (n N - 1 ), (4) 

7V=1 

+00 

f/T = J2 e^vMr^- 1 ), (5) 

N=l 
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Figure 3. Effective four- (a) and five-body recombination rates (b). The green dashed 
curve and the blue solid line follow the theoretical model for L4 and L 5 , respectively, 
with an additional scaling factor for L5; see text. The error bars include the statistical 
uncertainties from the fitting routine, the temperature and the trap frequencies. 



where represents the iV-body rate coefficients. The averaged atom densities are 
evaluated as (Z" 1 ) = fn N d 3 r = N N - 1 N-*l 2 [{mu 2 ) / (2iik B T)}^ N -V' 2 from a thermal 
distribution of the harmonically trapped atoms. Equation (jSJ) incorporates the anti- 
evaporation heating jl5] that results from the higher recombination rate in the densest 
part of the cloud (2e^ = 1 — 1/N). 

Under our experimental conditions, the first two terms of the sums can be omitted. 
One-body losses, resulting from collisions with the background gas, are negligible on our 
experimental time scale (Li = 0), and the use of atoms in their lowest sub-state assures 
that two-body losses vanish (L 2 = 0). Therefore, the three-body rate coefficient L 3 is 
the leading term contributing to the atom losses, and its resonant behavior is related 
to the Efimov trimers. This contribution is well understood and the coefficient L3 can 
be described by the well-established result of effective field theory jl6]. The rates of 
recombination events involving more particles are generally smaller than the one related 
to three-body losses and the contributions are difficult to discriminate because of the 
very similar behavior. Since the rate of recombination events for typical gas densities 
decreases rapidly with N, contributions with N > 5 are considered negligible in the 
following. 

A general fit to the experimental decay curves with L3, L4, and L5 as free parameters 
in Eq. (jl]) and Eq. ([5]) turns out to be practically impossible. Therefore, we fix L 3 
according to effective field theory, with parameters a_ = —955 ao and r\- = 0.08 as 
determined for three-body recombination (N = 3) in our previous experiment jl3]. We 
can now interpret the additional losses in terms of four-body and five-body decay. In 
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order to avoid any fitting ambiguities, we chose the simple approach to describe these 
losses either in terms of an effective four-body loss coefficient L 4 c fj (setting L 5 = 0) or 
an effective five-body loss coefficient L 5je fj (setting L 4 = 0). 

Figure |3]^a) shows L^g as extracted from our experimental data in comparison 
with the theoretical predictions, obtained by numerically evaluating L 4 from Eqs. (00) 
and fl3]) for our experimental conditions. Here we have adjusted the decay parameter 
to ?7_ = 0.33, which as a non-universal parameter depends on N. The comparison 
shows that the losses observed around —450 ao can be fully attributed to the four- 
body recombination resonance. The four-body loss peak position a 4j _ = —440(10) ao 
corresponds to 0.46(1) a_ and is in very good agreement with the theoretical value 
0.44(1) a_ [34] and previous observations [12]|TJ3[20l[39] . In contrast the enhancement of 
losses centered at about —300 ao cannot be interpreted in terms of the known universal 
four-body cluster stated^, suggesting that a different loss mechanism is present. 

An alternative representation of the same data in terms of L 5 e g is shown in Fig.|3]^b), 
together with the results of our theoretical model for L5; here we adjusted the relevant 
decay parameter to r/_ = 0.20. The model nicely explains the loss rates in the region 
where three- and four-body losses cannot account for the experimental observations. 
Remarkably, the resonance position = 0.64(2) a 4j _ is in agreement with the 

theoretical predictions 0.65(1) a 4] _ [36l[37]. However, quantitatively, the experimental 
values for L5 are about 15 times larger than the calculated ones. To account for this, we 
introduce a corresponding scaling factor. We find that this deviation might derive from 
non-universal effects that modify the value of the calculated WKB phase 7 by about 
10%, which remains within the realistic uncertainty range of our theory. 

An experimental search for higher-order recombination resonances (N > 5) that 
would be expected at lower values of the scattering length did not show clear signatures. 
A general problem arises, namely that the phase-space density cannot be further 
increased without causing the collapse of a Bose-Einstein condensate at negative values 
of a. To induce faster losses, adiabatic compression can increase the density, but then 
the higher temperatures cause increasing problems with the unitarity limit for high N. 
By decreasing the temperature, constraints by the unitarity limit can be avoided, but 
then losses for high N get so small that they become practically unobservable. 

Based on the above results and parameters, we model the general loss behavior 
in the region of interest. Figure |4] shows an example for three, four- and five-body 
recombination by plotting the atom losses, under typical experimental conditions, for 
a fixed hold time and variable scattering length. The "family portrait" of iV-body 
recombination highlights the different contributions and confirms how the different loss 
features dominate the losses at the resonant positions. The experimental data plotted 
in the inset show that the peak positions and the magnitude of losses are in very good 
agreement with the simulated losses. Note that the somewhat higher experimental 
losses can be attributed to an additional loss that occurs during the ramp to the target 

% We cannot rule out the possibility of a non-universal few-body state, which may be associated with 
higher partial waves. However, such an accidental coincidence appears to be rather unlikely. 
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Figure 4. Calculated and measured fraction of loss atoms from an atomic sample of 
initially 5 x 10 4 atoms at a temperature of 80 nK after a hold time of 100 ms. The 
red dotted line corresponds to the losses predicted for three-body recombination only, 
while the dashed green line and the blue solid line include also contributions from four- 
and five body recombination, as quantified in this work. A cut-off to the maximum 
losses has been applied according to Eq. $2$, as suggested in \$7\. The inset shows a 
close up of the region of dominant four- and five-body recombination and compares 
the theoretical results with the experimental observations. 



magnetic field strength. 
4. Conclusion 

By pushing the limits of ultracold few-body physics, we have explored a universal five- 
body recombination resonance both experimentally and theoretically. The observed 
series of recombination features, which we interpret in terms of three-, four- and five- 
body recombination resonances, provides crucial evidence for the existence of a family 
of universal iV-body bound states tied to Efimov trimers. The infinite series of iV-boson 
cluster states represents a paradigm for the general implications of Efimov physics for 
many-body systems. We speculate that similar scenarios also exist for other few-body 
systems of increasing size, containing fermionic constituents or particles of different 
masses, with important consequences for the interaction properties of the many-body 
system. 
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Appendix A. Extraction of the trapped five-body spectrum 

A crucial ingredient of our recombination analysis is the accurate extraction of the 
five-body trapped spectrum. The trapped few-boson Hamiltonian is given by 

N ( h 2 V 2 mcu 2 r 2 \ 
H = E + ^ +E^(r y -) + £ V 3b (R ijk ) (A.l) 

i=i \ zm z / i>j i>j>k 

where u is the trapping frequency, V 2 b{r) = V2boe~ r2 ^ 2r & is the two-body potential 
and V3b(R) = V3boe~ R2 ^ 2R ^ is a three-body potential. Here, r and Ro are the 
ranges of the two and three-body potentials that are fixed during the calculation, and 
Rijk = \J (rfj + rf k + r 2 k )/3. The two-body potential strength V2M is used to tune the 
two-body scattering length and V3& is set to zero in most calculations (see discussion 
below). To relate the model potential with the experimental we rely on the universal 
character of the low energy few-boson physics and we relate r with the van der Waals 
length r v d w so that the three- and four-body recombination peaks coincide with those 
experimentally observed. In this transformation, we relate the energy scales so that 
they scale as inverse length squared. 

To extract the recombination parameters, the five-body spectrum is analyzed as a 
function of the scattering length and the trapping confinement. The scattering length is 
changed in the region — 4tq < a < and the trapping frequency is changed so that the 



corresponding trapping length = yh/(muj) varies in the range 5r < a ho < 100r . 
This region of scattering lengths and energies corresponds to the low energy region where 
the five-body resonance occurs. 

The three-body interaction is used to shift the position of the five-body resonance 
and explore the dependence of the recombination on the scattering length. The range 
of the three-body potential is taken to be Rq, which is smaller than the two-body r so 
that it is mainly relevant at small hyperradii, i.e. so that its main contribution in our 
recombination formula is to change the inner phase. For the three-body interactions 
considered, the spectrum follows the linear dependence with the scattering length that 
is expected at small and negative scattering lengths. This linear dependence arises 
from only two-body physics which, in the hyperspherical framework, is described by the 
long-range behavior of the potential curves. Using the zero-range model of the two-body 
interaction, one can derive a first order correction of the hyperspherical potential curve. 
Here, we follow a similar procedure to that on Ref. [IE], but for a set of coordinates in 
which the center of mass has been removed. In this approximation, the lowest potential 
takes the form 

3(N - 2)(3iV - 4)h 2 y N u 2 R 2 h 2 a(N - 1)N^T (jN - |) 

V[R) ~ 8^ + + V2W23 r p _ 3) ' ( } 

The first two terms correspond, respectively, to the hyperangular or "mock-centrifugal" 
kinetic energy and the trapping potential; and the third term represents the interaction 
corrections. Equation IA.2I is only valid for small \a\ in the region where R is much larger 
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Figure Al. (a) The lowest two eigenenergies of a trapped five body system are shown 
as functions of the scattering length for different trapping frequencies. Different colors 
represent different trapping frequencies. The combination of these states essentially 
describes the energy of the five-body state in the inner region of the potential E mo i (a) 
(the diagonal curve). Here E sr = ft 2 /(mrp) and ro is the characteristic range of the 
two-body model potential that can be tuned to obtain the five-body resonance (i.e. 
ro ~ l-7r v d w where r v d w is the van der Waals length), (b) The near-threshold behavior 
of A. The fitting of the lowest energy points leads that A oc AE h . The lowest three 
points lead to b w 5.004 as expected from the known threshold behavior |40j . 



than the interaction range. Using perturbation analysis, one derives the well-known 
corrected energy of the lowest state, 

^ [2 a N(N - 1) /a \ 

E « 3 iV - 1 — + J tUj-^- '-. A.3 

2 V 7t a ho 2 

Our numerical calculations, with and without the three-body forces, show the linear 
behavior described in Eq. IA.3I in the region \a\ < |as,-|. This suggests that the 
hyperradial potential in the R ^> r is well described by Eq. IA.2I and that, in this 
region, the potential is independent of the three-body forces. Thus, it is consistent to 
interpret the main contribution of the three-body interaction as a modification of the 
short range physics that controls 0j n . 

The calculations are carried out using a correlated Gaussian basis set expansion 
limited to describe L p = + trapped states since the energetically lowest scattering 
continuum corresponds to zero angular momentum and positive parity. In this 
computation, we are only interested in the lowest trapped states which, in the 
hyperspherical picture, are supported mainly by the lowest potential curve. Therefore, 
the basis set is designed to accurately describe those states. The calculations include 
thousands of basis functions which are optimized for different scattering lengths and trap 
lengths. To verify the convergence of the energies, we carry out several optimization 
steps. The typical spectrum obtained by this analysis is shown in Fig. lAl( a). 
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Appendix B. Semiclassical analysis of the trapped spectrum 

To extract the semiclassical recombination parameters, we repeat the same semiclassical 
analysis used to derive Eqs. (2,3) but for a hyperspherical potential with the trapping 
potential (see Fig. |2]). This basically amounts to solve a double well problem using a 
semiclassical analysis. For this analysis, we follow closely the prescription in Ref. [12] 
to determine the quantization condition 

f3 = - tan -1 Q e ~ 27 tan(0 in )^ - 4> out + ^ = nn (B.l) 

where <pi n is the phase in the inner region, 7 is the barrier phase and <p ou t is the phase 
in the external trapped region (see Fig. [2]). The semiclassical phases are 

q(R)dR, (B.2) 

inner well 



7(£,a) = Im / q(R)dR, (B.3) 

./barrier 



■'out 



q(R)dR, (B.4) 
trap well 



where q(R) = yj2fi[E — V(R)\ and V(R) is the hyperspherical curve with the Langer 
correction. Here, the different integration regions are bounded by the classical turning 
points, i.e. the R positions at which q(R) = 0. In our analysis, <fi and 7 are assumed 
to be unaffected by the trapping potential which is expected to be negligible at small 
hyper radii. 

After some mathematical manipulation, the quantization condition can be written 

as 

-, — r- — 7 7 = - ( B -5) 

tan(0 in ) tan((/w) 2 

where A = e _27 /2. For tunneling energies well below the barrier height A 1 which 
implies that the quantization condition is fulfilled when either <fr in « vr(zj n + 1/2) or 
<Pout ~ ?r(w + 1/2). These conditions (<p in « n(i in + 1/2) and <p out « ir(i out + 1/2)) 
can be interpreted as having a bound state either in the inner or outer region. We call 
E in (a) and E out (a) the energies at which <p in = n(i in + 1/2) and 4> out = n(i out + 1/2), 
respectively. Here i a are integers representing the number of bound states supported by 
the inner and outer well respectively. Thus, close to those energies, the inner and outer 
phases take the form a = ir(i a + 1/2) + (p' a (E — E a (a)) where <f>' a = d<fi a /dE. In the 
vicinity of the avoided crossings, the quantization condition can be approximated to 

(E in (a) - E){E wt {a) - E) - A « 0. (B.6) 

Equation IB. 61 resembles a determinant of a 2 x 2 Hamiltonian matrix which lends 
itself to the following physical interpretation. The states supported by the inner and 
outer regions are coupled to each other and the coupling energy Vi nfiUt is related to the 
tunneling through the barrier by V^ nout = 2 A ,, — . The eigenenergies of this double well 
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problem (E\ and E 2 ) reach their squared minimal difference at the avoided crossing, 
equal to: 

(E 2 -E l ) 2 = 2A/(<p' m <p' out )- (B.7) 

To verify the validity of this equation, the energy dependence of A is explored by 
varying the strength of the trapping potential. According to the threshold laws |40j . 
for a five-boson problem, A increases with energy as A oc E 5 . Figure lATT b) confirms 
that the scaling of A calculated through the avoided crossing analysis shows excellent 
agreement with the threshold law prediction. 
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